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CALCULATED DIPOLE MOMENT AND ENERGY IN COLLISION OF A 
HYDROGEN MOLECULE AND A HYDROGEN ATOM 
by R. W. Patch 
Lewis Research Center 

SUMMARY 

Ab initio calculations of the interaction energy and electric dipole moment of -a' hy- 
drogen molecule (Hg) colliding with a hydrogen atom (H) were carried out on a digital 
computer using three Slater -type Is orbitals in the orthogonalized valence-bond theory 
of McWeeny. Each orbital exponent was optimized. The internuclear distance in the 

H„ molecule was varied from 7. 41599X10" 1 * to 7.67292x10“^ meter (1.401446 to 

* -10 
1.450000 bohrs). The intermolecular distance was varied from 0.529167x10 to 

2. 11667x10”*° meter (1 to 4 bohrs). Linear, scalene, and isosceles configurations were 
used. A weighted average of the interaction energies of the three equilibrium configura- 
tions was then taken for each intermolecular distance. 

The interaction energies obtained appear to be qualitatively correct but are higher 
than a number of published curves. This is because of the simple model used. The 
electric dipole moment should not be so sensitive to the model. 

The electric dipole moment and its derivative with respect to Hg internuclear dis- 
tance were consistent with classical values at large intermolecular distances with one 
minor exception, which is explained. Both sets of values are intended for use in calcu- 
lating pressure -induced vibrational absorption coefficients in the infrared. 


INTRODUCTION 

In high -temperature propulsion devices such as gas-core nuclear rockets, an im- 
portant mechanism of heat transfer is radiant energy exchange between volumes of gas 
and between the gas and the wall. When such devices use high-pressure hydrogen, it is 
necessary to know the strength of pressure -induced infrared absorption to calculate the 
heat transfer. This strength depends on the interaction energy and electric dipole mo- 
ment of the two colliding chemical species. The strength has been calculated for two 


colliding hydrogen molecules (H 2 ) (ref. 1) but not for H 2 colliding with a hydrogen atom 
(H). For this second case, which is the motivation of this report, the interaction energy 
is fairly well agreed upon (ref. 2), but the electric dipole moment had never been calcu- 
lated. For large intermolecular distances the electric dipole moment can be found from 
the quadrupole moment of H 2 and the polarizability of H, but for small intermolecular 
distances the Hg wave function must be known. 

The object of the work described herein was to calculate the electric dipole moment 
and its derivative and energy after determining the Hg wave function. The energy was 
desired only for comparison with experiments and other calculations. 

Past work on H 2 -H interaction may be divided into ab initio calculations, semi- 
eifapirical calculations, and experiments. These three types of work will be discussed 
in subsequent paragraphs. 

Ab initio calculations for H 2 -H or Hg have been carried out by a large number of in- 
vestigators (refs. 3 to 31). Much of the early work utilized rough approximations for 
some of the molecular integrals, making it of little contemporary value. This pitfall 
was avoided by Hirschfelder, Eyring, and Rosen (ref. 3), Hirschfelder , Diamond, and 
Eyring (ref. 4), and Hirschfelder (ref. 6), who exhausted a basis of three Is atomic 
orbitals centered on the protons. More recently, Shavitt, Stevens, Minn, and Karplus 
(ref. 27) exhausted a basis of 15 orbitals (Is, Is’, 2p x , 2p y , 2p z on each proton). This 
work was extended by Porter, Stevens, and Karplus (ref. 30) . The most accurate work, 
however, is probably that of Conroy and Bruner (ref. 26), who improved the electron 
correlation by using a wave function that was the product of a shape series, a singlet- 
pair correlation series, and a triplet-pair correlation series. 

A large number of semiempirical calculations for H 2 -H or Hg have also been made. 
For large inter molecular distances the treatments of Margenau (ref. 32) and Mason and 
Hirschfelder (ref. 33) are noteworthy. For small intermolecular distances, the treat- 
ment of Porter and Karplus (ref. 2) correlates experimental data best. None of these 
are useful in finding electric dipole moment. 

Experimental studies of H 2 ~H interaction are also numerous. Amdur (ref. 34) 
measured the scattering of an H-atom beam in H 2 gas. The potential he obtained is 
.iower than that of Porter and Karplus (ref. 2). Schulz and LeRoy (ref. 35) measured 
' the rate for the reaction 


H + p-H 2 - o-H 2 + H 

in a flowing. system. Karplus and Porter (ref. 36) compared this with the potential of 
Porter and Karplus (ref. 2) and found fair agreement. Geddes, Krause, and Fite 
(ref. 37) used modulated crossed beams to obtain the reaction 
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D + H 2 -HD + H 

They measured velocities and angular distribution of HD. Brumer and Karplus (ref. 38) 
found that the results of reference 37 agreed with the potential of Porter and Karplus 
(ref. 2). 

In this study no attempt was made to do the most sophisticated possible ab initio 
calculation for H 2 - H because accuracy within a few percent for the dipole moment was 
all that was desired. To achieve this, a full valence-bond ab initio calculation for 
H 2 - H was carried out with accurate values for all molecular integrals. In order to 
cover nuclear configuration space, linear, scalene, and isosceles configurations (see 
fig. 1) were included with inter molecular distances from 0. 529167X10 -10 to 
2. 11667X10 -10 meter (1 to 4 bohrs). In three of the configurations (figs. 1(a), (c), 
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(a) Linear configuration A with equilibrium internuclear distance. 
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(b) Linear configuration B with greater-than-equilibrium inter- 
nuclear distance. 



(c) Scalene configuration C with equilibrium internuclear distance. 



(d) Scalene configuration D with greater-than-equilibrium inter- 
nuclear distance. 
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(e) Isosceles configuration E with equilibrium internuclear distance. 
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(f) Isosceles configuration F with greater-than-equilibrium inter- 
nuclear distance. 


Figure 1. - Six configurations of a hydrogen molecule and a hydrogen 
atom used in calculations. I ntermolecular distance R r had values 
from 0. 53x10“*® to 2. 12x10“*® m (1 to 4 bohr). I nternuclear dis- 
tances in the figure are in bohrs (1 bohr = 5. 29167x10“** m). 
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and (e)) the H 2 molecule had equilibrium internuclear distance. The average H 2 - H 
interaction energy was then calculated for a given inter molecular distance by taking 
weighted averages of the interaction energies of these three configurations. In the other 
three configurations (figs. 1(b), (d), and (f)) the H 2 molecule had a greater -than- 
equilibrium internuclear distance so that the derivative of the dipole moment with re- 
spect to internuclear distance could be calculated from the electric dipole moment of all 
six configurations. 

Professor V. Magnasco of the Istituto di Chimica Industrial delTUniversita, 

Genoa, provided listings and a deck of various computer programs for evaluating molec- 
ular integrals. Although these programs were not used for this report, they were use- 
ful for checking and provided ideas for the programs that were used. 

ANALYSIS 

Model and Configurations 

The interaction of H 2 and H may be. calculated by the orthogonalized valence-bond 
theory of McWeeny (refs. 39 and 40), which utilizes symmetrically orthogonalized 
atomic orbitals in the Born-Oppenheimer approximation. In this method orthogonalized 
orbitals are formed from linear combinations of atomic orbitals. These are multiplied 
by spin eigenfunctions to give orthogonalized spin orbitals. Formal orthogonalized 
valence -bond structures are formed from linear combinations of anti symmetrized 
products of the orthogonalized spin orbitals. The system wave function consists of 
linear combinations of these structures. This involved procedure merely amounts to 
taking the proper linear combination of the hextuple products of three atomic orbitals 
and three spin eigenfunctions. This method allows full configuration interaction without 
unduly complicating a digital computer program and allows straightforward calculation 
of the electric dipole moment. It has previously been successfully applied to the inter- 
action of two Hg molecules by Magnasco and Musso (refs. 41 and 42) and Patch (ref. 43). 

The model chosen had a basis of three Is Slater-type atomic orbitals Xp centered 
on the protons and had full configuration interaction. The three normalized Is atomic 
orbitals were 

AsV / 2 

Xp(i) =[-£j e p pl p = a,b,c s (1) 
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where the i in parentheses indicates the i^ electron, and r . is the distance be- 
tween proton p and electron i (symbols are given in the appendix). An orbital 
Xp(i), of course, reduces to the wave function for the ground state of an H atom if 
? = 1,88976x10^ reciprocal meters (1.0 bohr"*). 

If all the atomic orbitals were assigned the same orbital exponent C, the model 
would not accurately represent H 2 + H at large intermolecular distances because isolated 
Hg and H have different optimum orbital exponents. If the two orbital exponents belong- 
ing to orbitals centered on H 2 protons were assigned the value for isolated H 2 , and the * 
other was assigned the value for isolated H, inconsistencies at small intermolecular dis.- 
tances would result (the linear symmetric configuration would have different orbital ex- 
ponents for the end protons, and the equilateral triangle configuration would have unequal 
orbital exponents). Hence, each orbital exponent was optimized for each intermolecular 
distance and each H 2 internuclear distance of each configuration. 

In ab initio calculations the 'average interaction energy should be based on the energy 
of the complex minus the energies of the two separated chemical species (in this case H 2 
and H) calculated in a consistent manner. In the limit of infinite intermolecular dis- 
tance, the model in this report results in an H atom and a molecule (H 2 ) described by the 
covalent -ionic valence-bond model of Weinbaum (ref. 44). Because of the ultimate ap- 
plication of this report to pressure -induced absorption, the experimental equilibrium 

internuclear distance of 7. 41 599x10“ 11 meter (1.401446 bohrs) for H„ (ref. 45) was used 

-11 ^ 

rather than Weinbaum's value of 7.4962x10 meter (1.4166 bohrs). Hence an isolated 

1 8 

H„ molecule had an energy E of -5.0037x10 joule (-1 . 14779 hartrees) rather than 
Weinbaum's value of -5.0042x10“ joule (-1.14790 hartrees) (ref. 43). The average 
interaction energy is found by averaging over nuclear configurations, which was facili- 
tated by choosing the configurations shown in figures 1(a), (c), and (e). The other three 
configurations possessed a greater -than-equilibrium internuclear distance. 

In all six configurations the intermolecular distance R r was varied. It had values 
of 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.8, 2.0, 2.5, 3.0, 3.5, and 4.0 bohrs (1 bohr = 
5.29167X10' 11 m). 


v Energy and Interaction Energy 

To find the electronic energy of the H 2 - H complex, a Lowdin transformation was 
performed on the orbitals, resulting in orthogonalized orbitals. Formal orthogonalized 
valence-bond structures were formed, and the matrix elements of the Hamiltonian were 
found. An eigenvalue problem was then solved to obtain the electronic energy. The nu- 
clear repulsion energy was added to the electronic energy to get the energy of the 
complex. The energies of isolated H and H 2 were subtracted from the energy of the 
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complex to yield the interaction energy. This entire procedure was repeated for differ- 
ent values for the orbital exponents to enable the three orbital exponents to be optimized 
by numerical methods. These steps are given in the following subsections. 

Lowdin transformation . - The symmetric orthogonalization procedure of Lowdin 
(ref. 46) was applied to the y • The overlap integral is 


S pq*/V 1 N <1)dv l 


( 2 ) 


Let 


w 

pq 2 w 


(3) 


where R is the distance between protons p and q. If fL = (! q , S is given by 
Slater (ref. 47). If C p ^ ? q , 



e q -e“ q 2 

pq 4 

T w ^ 

L pq 


w. 


w + 

pq 


H| + e- ( MU-2- + -2_ 

q q 


“ i 2 i 3 


-W 


pq 


w. 


pq 


(4) 


where 


C - If 

q = _P _SL R 

2 ^ 


(5) 


It is desired to find the Lowdin transformation matrix such that 

X. 


xi= x M ^ 1/2 


( 6 ) 


where x is a row vector with elements Xp> and xl is a row vector of the three 
orthogonalized orbitals First the eigenvalue problem 
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(7) 


USU = d 


is solved, where S is a matrix with elements S . The matrix d is diagonal, with 

pq 

eigenvalues as diagonal elements. Thus 


,- 1/2 _ 


,- 1/2 

'll 

0 


0 

j— 1/2 

*22 

0 


0 

0 

,- 1/2 

'33 


( 8 ) 


Finally, 


M“ 1//2 = Ud“ 1//2 U 
X 


(9) 


Use of equation (6) then accomplishes the Lowdin transformation to obtain x\ 

Formal orthogonalized valence -bond structures . - To produce the proper anti- 
symmetry, the are combined in antisymmetrized products of orthogonalized spin 
orbitals. 


(pqF) = (3’.r 1/2 det[x^(l)a(l),x^(2)a(2),x;(3)|3(3)] 


( 10 ) 


where a(l) is the spin eigenfunction of electron 1 with the component of spin angular 
momentum along the axis of quantization equal to R/2, and /3(3) is the spin eigenfunction 
of electron 3 with the component of spin angular momentum along the axis of quantiza- 
tion equal to -fi/2. The shorthand notation is given on the left side of equation (10), 
wherein a bar indicates /3 spin and no bar indicates a spin. To facilitate computing 
for some purposes, the (pqr) are arranged in dictionary order in a three by three matrix 
(ref. 48): 


(aba) 

(abb) 

(abc) 


(aca) 

(acb) 

(acc) 

(ID 

(bca) 

(bcb) 

(bcc) 



For our purposes we will regard the (pqr) as a nine -element row vector with the ele- 
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ments numbered down the first column in equation (11), then down the second column, 
etc . When numbered in this way , set 

<p i = (pqr) i = 1, 2, . . . , 9 . (12) 


where p, q, and r are given in equation (n) . 

The formal orthogonalized valence-bond structures p^ are linear combinations of 
the (p^. These p^ are eigenfunctions of the total electron spin with total spin quantum 
number of 1/2 and component of the total spin angular momentum along the axis of quan- 
tization of R/2. These are obtained by temporarily putting a phantom orbit at in- 
finity as outlined by Pauling (ref. 49) and are 


1 ~ 2 1//2 (-<P 3 - cp§) 

^ 5 = <?1 

> 2 = 2~ 1/2 (-<P 3 + < Prj) 

05 

II 

1 

CO 

'3 = ~^2 

^7 = <P 4 

'4 = “*6 

P Q =<P 8 





(13) 


where formal singly polar structures are included. All jf/. are orthonormal (see 
ref. 50) except p^ and p 2 • These may be made orthonormal by the Schmidt ortho- 
gonalization procedure, which results in replacing p 2 i n equation (13) by a new p 2 
such that 


p 2 = 6" l//2 (-<p 3 + <p 5 + 2c prj) (14) 

Equations (13) (excluding p 2 ) an< ^ (14) can be represented by 


p = <pV_ (15) 

where V is a nine by eight matrix. The wave function of the system was taken as 

* = P £ (16) 

where C is a column vector determined by minimizing the electronic energy E e _ This 
is accomplished by the variation principle for orthonormai functions (ref. 51), which 
poses an eigenvalue problem. 
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( 17 ) 


HC = E C 
— e — 

where the smallest of the several eigenvalues E g is the ground state electronic energy 
and H is a matrix whose elements are given in the following sections. 

Matrix elements of the Hamiltonian . - The Hamiltonian operator for the electronic 
energy in atomic units with the hartree as the unit of energy is 


i-1 


H = 


x-^XX-XXXt 

1_i i=l p=a i=2 j=l 


(18) 


2 

Here is the Laplacian operator for the coordinates of electron i, and the p sum- 
mation is over protons a, b, arid c. The elements of the H matrix in equation (17) 
are then 


V'/ *l**k a 7 , < 19 ’ 

where r includes electron configuration and spin coordinates. From equations (15) 
and (19) 

H = VH^V (20) 

where 


( h 4j = / 'V'-'j dT 1211 

The matrix elements IH ) can be found from Slater's rules (ref. 50 after correction of 
a typographical error) if two kinds of integrals are known. 


X 

P =a 
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(23) 


J ijkZ 


/ j X , i (l)Xj(l)X , k (2)X;(2) ~ 


dvj dvg 


where is the distance from electron 1 to electron 2, dv^ is the volume element for 
electron 1, and dvg is the volume element for electron 2. From equations (6), (22), 
and (23) 


G ijkZ 


where 


- £ « /2 ) ki £ k* ♦ X 

k=l 1 = 1 \ P- a 

Z k i/ 2 lZ ( m * 1/ % Z k / 2 lZ m 


p z Gmn °P 


m=l 


n=l 


o=l 


p=l 


K k I’f 

L ikp m J X i <1)Xk<1) (i“Z 


mnop 


if 


x m (1)x n (1)x o (2)x p (2) ~ dv l dv 2 

Sr r» 
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(24) 


(25) 


(26) 


(27) 


(28) 


The methods of evaluating , L j k p, and G mno p are given by Patch (ref. 43). 
Energies. - The energy E is the electronic energy plus the nuclear repulsion 


energy : 


E = E Q + _L+_L_ + _L 

e 


^ab ^ac ^bc 


(29) 
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The interaction energy E^ ^ is the energy E minus the energies E m (r) and E^ ' of 
the isolated 
spectively: 


the isolated Hg molecule with internuclear distance r and the isolated H atom, re- 


E int = E - E m< r > - E at < 30 > 

The diatomic internuclear distance used in equation (30) was the same as in the Hg - H 
complex (see fig. 1), and the corresponding E m (r) have been tabulated by Patch 
(ref. 43) for Weinbaum’s covalent -ionic valence-bond model of Hg. Of course, E a ^ 
was -2 . 17972X10" 18 joule (-0 . 5 hartree) . 

Optimization of orbital exponents . - To optimize the orbital exponents, it is neces- 
sary to solve the eigenvalue problem (eq. (17)) for a number of sets of orbital expo- 
nents. A least-squares fit is made to the resulting values of E^, and the values of the 
orbital exponents that minimize E^ are found from the fit coefficients. 

The simplest case involves isosceles configurations because by symmetry £ b = £ 
(see fig. 1(d)). If the optimum values ? aQ and £ bo are known approximately, then a 
least-squares fit may be made. 


E int 


W o + W !<a + w 2*b + W U«a + W 22«b + W l 2 Vb 


(31) 


where the E int to be fit are calculated from equations (17), (29), and (30) for values 

?a - ?ao - ?ao' ?ao * A ?. and ?b = ?bo ' A ?- W ?bo + A ? < nine sets ot ?a 
and ^ b ). Here, is an arbitrary number much smaller than £ ao or ^ Q , and £ aQ 

and ^ bo are estimates. After the least-squares fit coefficients W have been found , 

the minimum of equation (31) can be found by direct calculation, yielding the true £ „ 

dU 

and C bQ . 

The other case applies to linear and scalene configurations. Here, the optimum or- 
bital exponents are generally all different so the least-squares fit is of the form 


E int “ W o + W l?a + W 2?b + W 3?c + W ll<a + W 22^b + W 33«c + W l 2 Vb 


+ W 13 Vc + W 28 Vc 


(32) 


where the E.^ to be fit are calculated from equations (17), (29), and (30) for values 

%L - ? ao - ?ao> <,0 ♦ A ?; ?b - ?bo ' A? ' W ?bo * ?c A f co - <co> 

(27 sets of ^ a , ^ b , and £ ). The rest of the procedure is the same. 
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Average Interaction Energy 


The average interaction energy of an Hg molecule and ah H atom is, to a first 

approximation, the interaction energy for equilibrium diatomic internuclear distance av- 
eraged over all possible orientations of the Hg molecule. 




(33) 


where £2 is the solid angle for Hg orientation. 

Since only three configurations with equilibrium diatomic internuclear distances 
were considered in this report, an additional approximation was necessary. This was 
the truncated expansion of in terms of Legendre polynomials O^: 

E int = f O 0 O (cos ^ + f 2 0 2^ co ® ^ + f 4 ° 4 ( cos ^ (34) 


where the f’s are expansion coefficients and 6 is the angle between the internuclear 
axis of Hg and the line through the H atom bisecting Hg. Odd Legendre polynomials 
were ruled out by symmetry. The result of substituting equation (34) into (33) is 




int,A 


. 8 -r, , 2^ p 

l g int , C g int , E 


(35) 


Electric Dipole Moment 

Calculation of the dipole moment requires an electron population analysis in terms 
of the original nonorthogonal orbitals Xp* A column vector CV may be defined by 

* = (pV (36) 

From equations (15), (16), and (36) 


= VC . (37) 

The electron population in the orbital density X' m X' n irrespective of spin state is 
(ref. 41) 
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P- = V C’ V C’ (n mano1 + n m ^) (38) 

mn / v r / y s \ rs rs ) v ' 

r=l s=l 


where n ™ Q ' ncv and n™f n ^ are coefficients of spin-orbital products and are called n^ 
r s r s rs 

by Magnasco and Musso (refs. 41 and 48), who also give rules for evaluating them 
(ref. 48). 

The electron population in the orbital density X p X q irrespective of spin state is 
(ref. 41) 


pq 


ZK 1 / 2 lmZ^k 1/2 t 
m=l n=l 


(39) 


The normalized 1 -electron population is 



= P S 

pq pq 


(40) 


where the three diagonal elements of the matrix give the atomic populations and the 
off-diagonal elements give the overlap populations. For instance, the population in the 
overlap between orbitals centered on protons a and b is ^ a k + ^ba‘ The elements 
of the matrix, of course, add to three, the number of electrons. 

To find the dipole moment from ^ it is necessary to know locations of the atomic 
charge centers and the overlap charge centers. Cartesian coordinates with the origin 
one-third of the way from the molecule toward the atom were chosen (see fig. 1(a)). The 

i-L 

j Cartesian coordinate of either type of charge center was designated C^., where p 
and q correspond to the two orbitals x p and x q involved in an overlap charge or 
where p and q are equal and correspond to the one orbital involved in an atomic 
charge. The C pp j for atomic charges are, of course, identical to the coordinates of 
the protons C p „ For overlap charges 




Pi 



^pq) c qj 


(41) 


V 


where A is the displacement of the overlap charge center from the p - q geometric 
PQ 

center, measured in the direction from p towards q and expressed in units of R . 

K S p = ? q , A pq - 0 by symmetry. If C p * C q , 
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R;Le 

a = pq 
pq 


3 „- W pq 


8S pq w pq 




>-q -i- 


i i 


:2 , 


- e qfi.-L 

-2 


— + — + 3 f w. 


2 w 

lW pq Pq 


pq 


4 + J. + I + l 


q 2 q 3 q 4 


q 2 q 3 q 4 /J 




(42) 


IL 

The j component of the dipole moment in atomic units is then 

3 3 


v-ZZ 


5 C 

pq pqj 


p=l q=l 


(43) 


By symmetry, only the x and z components may be nonzero. 


Computer Program 

All calculations including optimizing the orbital exponents were done by one 
FORTRAN IV IBM version-13 program for an IBM 7094 digital computer. To obtain a 
set of optimized orbital exponents, interaction energy , and dipole moment, the running 
times were 4.5, 8.7, and 2.1 minutes for the linear, scalene, and isosceles configura- 
tions, respectively. 


RESULTS AND DISCUSSION 

In this section values for the orbital exponents, interaction energy, average inter- 
action energy, dipole moment, and dipole moment derivative are given and compared 
with the results of other investigators. 


Orbital Exponents and Interaction Energy 

In optimizing the orbital exponents, it was found that the estimates had to be fairly 
accurate for the coefficients of equation (31) or (32) to give meaningful values for the op- 
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TABLE II. - AVERAGE INTERACTION ENERGY 


FOR H 2 + H 


Intermolecular distance, 
R r 

Average interaction energy, 
R int 

m 

bohr 

J 

hartree 

5.29167X10 -11 

1.0 

2.5712X10 -18 

0.58980 

5.82084 

1.1 

2.0969 

.48100 

6.35000 

1.2 

1.7730 

.40670 

6.87917 

1.3 

1.4969 

.34336 

7.40834 

1.4 

1.2400 

.28445 

7.93751 

1.5 

1.0374 

.23797 

8.46667 

1.6 

0.8753 

.20079 

9.52501 

1.8 

.6367 

.14605 

1.05833X10' 10 

2.0 

.4732 

.10854 

1.32292 

2.5 

.2353 

.05398 

1.58750 

3.0 

. 1173 

.02690 

1.85208 

3.5 

.0569 

.01305 

2.11667 

4.0 

.0267 

.00613 




Interaction energy, E inl , hartrees 



Intermolecular distance, R r , bohr 
(a) Linear configuration A (see fig. l(al). 


10° p 


10 ' 3 ' — 



I ntermolecular distance, R r , bohr 


(b) Scalene configuration C (see fig. 1(c)). (c) Isosceles configuration E (see fig. 1(e)). 

Figure Z - Interaction energy for three configurations of a hydrogen molecule and a hydrogen atom, each configuration having equilibrium inter- 
nuclear distance in the molecule. 
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timum orbital exponents. Consequently, before running the program in the mode for 
automatic least-squares fit to and optimization of it was run in another mode 

that allowed one to be varied at a time, and the resulting E- nt 's were hand plotted. 
In this way good estimates for the were eventually obtained, allowing the program to 
be run in the automatic -fit mode. In some cases the automatic -fit mode had to be re- 
peated with progressively smaller A£ until a A£ of 1.89x10 reciprocal meters (0.01 
bohr~^) ran successfully, producing optimum t^'s that fell within the grid used to 
obtain the least-squares fit for that run. The final orbital exponents and interaction 
energy are given in table I for 3 nuclear configurations and 13 intermolecular distances. 
The apparent discontinuity at R r = 6.422x10"^ meter (1.2136 bohrs) for orbital expo- 
nents of configuration E (isosceles) is due to a change of symmetry of the eigenfunction 
as the isosceles triangle goes through the equilateral triangle configuration. 

The interaction energies are plotted in figures 2(a), (b), and (c) for configurations 
A, C, and E. The change of symmetry mentioned produces a slight Jahn-Teller cusp in 
the interaction energy of configuration E. For all configurations, the interaction energy 
increases monotonically as R r is decreased. For large R r the order of increasing 
interaction energy is configurations E, C, and A. For small R r the linear configura- 
tion A has the highest interaction energy. For intermediate R r it is difficult to gen- 
eralize. 


Average Interaction Energy 

The average interaction energy is given in table II and plotted in figure 3. It in- 
creases monotonically as R r decreases. 


Electric Dipole Moment and Its Derivative 

The components of the electric dipole moments for configurations A , C, and E are 
given in table in for 13 intermolecular distances. They are plotted in figure 4. Again, 
there is an apparent discontinuity due to the change of symmetry as the isosceles con- 
figuration goes through the equilateral triangle configuration. The x-component of the 
dipole moment for configurations A and E is zero by symmetry. 

To obtain the derivatives of the components of the electric dipole moment with re- 
spect to the Hg internuclear distance, all calculations were repeated for a greater -than- 
equilibrium internuclear distance (configurations B, D, and F of fig. 1). The derivatives 
were found from the approximation d/x^/dr « A/ju/Ar. The results are given in table III 
and plotted in figure 5. There is the usual apparent discontinuity for the isosceles con- 
figuration. 
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TABLE I. - ENERGY AND ORBITAL EXPONENTS FOR H 2 + H 


Intermolecular distance, 
R r 

Con- 

figura- 

tion 2 

Orbital exponent, 

Orbital exponent, 
« b 

Orbital exponent, 

Interaction energy, 
E int 

m 

bohr 

1/m 

1/bohr 

1 /m 


1/m 

1/bohr 

J 

hartree 

5.29167X10"** 

1.0 

A 

2.4966X10 10 

1.32113 

-3.7651X10 10 

1.99239 

1.8262x10*° 

0.96636 

8.7946X10"* 8 

2.01737 



C 

.2.4421 

1.29230 

3.1230 

1.65260 

1.8172 

.96159 

2.3052 

0. 52878 



E 

3.1118 

1.64667 

1.9325 

1.02262 

1.9325 

1.02262 

1.8885 

.43321 

5.82084X10'** 

1.1 

A 

2.4870X10 10 

1.31605 

3.6068x10*° 


1.8652x10*° 

0.98700 

5.5201X10'* 8 

1.26625 



C 

2.3651 

1.25154 

3.0987 

1.63971 

1.8481 

.97796 

1.8935 

0.43434 



E 

3.0355 

1.60630 

1.9094 

1.01037 

1.9094 

1.01037 

1 . 7975 

.41233 

6.35000X10' 11 

■ m 

A 

2.4567X10 10 

1.30003 

3. 4497X1 0 10 

1.82548 

1.8951x10*° 

1.00284 

3.6986X10"* 8 

0.84841 



C 

2.2972 

1.21560 

3.0642 

1.62145 

1.8765 

0.99298 

1.5612 

.35811 


1 

E 

2.9606 

1.56663 

1.8866 „ 

0. 99831 

1.8866 

.99831 

1.7344 

.39786 

6.87917X10"** 

■ 

A 

2. 4145X10 10 

1.27769 

3. 31 39X1 0 1 0 

1.75359 

1.9199X10*° 

1.01596 

2.5893X10"* 8 

0. 59395 



C 

2.2364 

1.18341 

3.0197 

1.59791 

1.9039 

1.00749 

1.2967 

.29744 


1 

E 

1.7928 

0. 94870 

2.3230 

1.22928 

2.3230 

1.22928 

1.5816 

.36281 

7.40834X10' 11 

1.4 

A' 

2.3658X10 10 

1.25189 

3.1988X10 10 

1.69272 

1.9416X10*° 

1.02742 

1.8747X10'* 8 

0.43003 



C 

2.1818 

1.15454 

2.9724 

1.57288 

1.9294 

1.02100 

.1.0868 

.24930 



E 

1.8188 

0. 96243 

2.3199 

1.22763 

2.3199 

1.22763 

1.3386 

. 30705 

7.93751X10' 11 


A- 

2. 3146X10*° 

1.22481 

3.1002X10 10 

1.64052 

1.9612x10*° 

1.03780 

1.3964X10"* 8 

0.32032 



c 

2.1334 

1.12894 

2.9216 

1.54600 

1.9537 

1.03382 

9.1962X10'* 9 

.21095 



E 

1.8344 

0. 97072 

2.3177 

1.22647 

2.3177 

1.22647 

1.1346X10'* 8 

.26026 

8.46667X10' 11 

1.6 

A 

2.2639X10 10 

1.19799 

3.0136X10*° 

1 . 59469 

1.9794X10*° 

1.04744 

1.0672X10"* 8 

0.24481 



C 

2.0912 

1.10657 

2.8681 

1.51773 

1.9768 

1.04607 

7.8522X10"* 9 

.18012 



E 

1.8444 

0.97600 

2.3155 

1.22530 

2.3155 

1.22530 

9.6348 

.22101 

9.52501x10"** 

1.8 

A 

2. 1708X10 10 

1. 14871 

2.8651X10*° 

1.51611 

2.0133x10*° 

1.06535 

6.6961X10"* 9 

9 



C 

2.0241 

1.07111 

2.7571 

1.45896 

2.0210 

1.06944 

5.8578 




E 

1.8582 

0. 98332 

2.3093 

1.22200 

2.3093 

1 . 22200 

6.9912 


1.05833x10"*° 


A 

2.0940X10 10 

1.10808 

2.7389x10*° 

1.44931 

2.0455X10*° 

1.08239 

4.5787X10"* 9 

0. 10503 



C 

1.9782 

1.04678 

2.6475 

T. 40097 

2.0638 

1.09209 

4.4641 

. 10240 



E 

1.8673 

0.98813 

2.3020 

1.21812 

2.3020 

1.21812 

5.1136 

. 11730 

1.32292x10'*° 


A 

1.9774X10 10 

1.04637 

2.4945x10*° 

1.32001 

2.1268x10*° 

1.12541 

2. 2612x1 (T 19 

0.05187 



C 

1.9252 

1.01876 

2.4308 

1.28630 

2.1625 

1.14432 

2.3275 

.05339 



E 

1.8844 

0.99714 

2.2849 

1.20908 

2.2849 

1 . 20908 

2.4029 

.05512 

1.58750x10"*° 

3.0 

A 

1.9314X10 10 

1.02201 

2.3544X10*° 

1.24587 

2.1983x10*° 

1.16326 

1.2 180X10."* 9 

0.02794 



C 

1.9101 

1.01077 

2.3260 

1.23083 

2.2240 

1.17686 

1.1875 

.02724 



E 

1.8942 

1.00237 

2.2749 

1.20382 

2.2749 

1 . 20382 

1.1452 

.02627 

1 . 85208x10"*° 


A 

1.9125x10*° 

1.01204 

2. 2973X10*° 

1.21567 

2.2382X10*° 

1.18440 

ISgBSBBl 

0.01432 



C 

1.9037 

1.00737 

2.2884 

1.21093 

2.2500 

1.19065 

lijSgi 

.01337 



E 

1.8970 

1.00381 

2.2706 

1.20155 

2.2706 

1.20155 


.01240 

2.11667x10'*° 

4.0 

A 

1 . 9029x10*° 

1.00696 

2.2781X10*° 

1 . 20548 

2.2554x10*° 

1.19347 

3.0254X10" 20 

0.00694 



C 

1 . 8990 

1.00491 

2.2756 

1.20415 

2.2604 

1.19613 

2.7552 

.00632 



E 

1.8960 

1.00328 

2.2690 

1.20070 

2.2690 

1 . 20070 

2.5023 

.00574 


a See fig. 1. 
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TABLE m. - COMPONENTS OF DIPOLE MOMENT AND THEIR DERIVATIVES FOR Hg + H 


Intermolecular distance, 
B r 

Con- 

figura- 

tion 3 

Component of dipole moment, 

Component of dipole moment, 
"z 

Derivative of component 
of dipole moment, 
3(i x /3r 

Derivative of component 
of dipole moment, 

m 

bohr 

C-m 

-at. units 

C-m 

at. units 

C 

at. units 

c 

at. units 



A 

0 

0 

2 . 67638 X 10 -30 

0.315693 

0 

0 

- 2 . 595 X 10" 20 




C 

1 . 81431 X 10' 30 


2.80947 

1 331392 

. 823 X 10' 20 

.0514 

- 1.554 




E 

0 

0 

2.12124 

.250212 

0 

0 

- 1.352 

-.0844 

5 . 82084 X 10" 11 

1.1 

A 

0 

0 

2 . 34716 XKT 30 

0.276860 

0 

0 

- 2 . 719 X 10' 20 

- 0.1697 



C 

1 . 69238 X 10" 30 

.199625 

2.52880 

.298285 

1 . 203 X 10" 20 


0.003 

. 9 



E 

0 

0 

2.48750 

.293414 

0 


- 2.022 


6 . 35000 X 10" 11 

1.2 

A 

0 

0 

2 . 09418 X 10' 30 


0 

0 


- 0.1319 



C 

1 . 57656 X 10' 30 

.185964 

2.17885 

.257007 

1 . 405 X 10" 20 

.0877 





E 

0 

0 

2.85553 


0 

0 


IBB 21 

6 . 87917 X 10' 11 


A 

0 

0 

1 . 85676 X 10 -30 

0.219015 

0 

0 

- 1 . 032 X 10' 20 

- 0.0644 



C 

1 . 45789 X 10' 30 

.171966 

1 . 77748 

.209663 

1 . 762 X 10 -20 


, 3.571 

.2229 



E 

0 

0 

- 4.18773 

r . 493965 

0 

0 

1.956 

.1221 

7 . 40834 X 10' 11 

1.4 

A 

0 

0 

1 . 611 16 X 10 -30 

0,190045 

0 

mm 

0 . 178 X 10' 20 

0.0111 



. C 

1 . 34361 X 10" 30 

.158486 

1 . 34934 

.159162 

1 . 964 X 10* 20 


5.120 

.3196 



E 

o 

0 

- 3.95336 

-.466320 

0 

del 

1 . 182 

.0738 . 



A 

0 

0 . 

1 . 35294 X 10" 30 

0.159586 

0 

0 • 

1 . 354 X 10' 20 




C 

1 . 23119 X 10' 30 

.145225 

9 . 18763 X 10" 31 

.108373 

2 . 107 X 10" 20 

.1315 

6.263 

.3909 



E 

0 

0 

- 3 . 75114 X 10' 30 

-.442467 

0 

0 

0.644 

.0402 

8 . 46667 X 10' 11 

1.6 

A 

0 

0 

1 . 08738 X 1 O' 30 

0.128262 

0 

0 

2 . 378 X 10' 20 

0.1484 



C 

1 . 12113 X 10' 30 

.132243 

5 . 08998 X 10' 31 

.060039 

2 . 185 X 10' 20 

.1364 

6.919 

.4319 



E 

0 

0 

- 3 . 55798 x 1 0 -30 

-.419683 

0 

0 

0 . 364 

.0227 

9 . 52501 X 10" 11 


A 

0 


5 . 74480 X 10' 31 

0.067763 

0 

0 

3 . 685 X 10" 20 

0.2300 



C 

9 . U 158 X 10" 31 


- 1.76507 

-.020820 

z.uaxio ' 20 

.1341 

6.825 

.4260 



E 

0 


- 3 . 14938 X 10' 30 

-.371486 

0 

0 

- 0.328 

-.0205 

1 . 05833 X 10" 10 

2.0 

A 

0 

0 

1 . 56746 X 10' 31 

0.018489 

0 

0 

3 . 872 X 10' 20 

0.2417 



C 

7 . 19679 X 10" 31 

.084890 

- 6.21218 

-.073276 

1 . 913 X 10" 20 

.1194 

5.466 

.3412 



E 

0 

0 

- 2 . 71693 X 10' 30 

-.320476 

0 

0 

- 0.692 

-.0432 



A 

0 

0 

- 2 . 47144 X 10" 31 

- 0.029152 

0 

0 

2.075 

0.1295 



C 

3 . 47420 X 10" 31 


- 8.27457 

-.097603 

1 . 043 X 10" 20 

.0651 

1.804 

.1126 



E 

0 

0 

- 1 . 66838 X 10' 30 

-.196794 

0 

0 

- 0.838 

-.0523 

1 . 58750 X 10" 10 


A 

0 

0 

- 1 . 59772 X 10' 31 

- 0.018846 

0 

0 

1.237 

0.0772 



C 

1 . 46979 X 10 -31 

.017337 

- 5.31294 

-.062669 

. 418 X 10 -23 

.0261 


.0334 


in 

E 

0 

0 

- 9.12345 

-.107616 

0 

0 


-.0341 

' 1 . 85208 X 10" 10 

H 

A 

0 

0 

- 6 . 15572 X 10' 32 

- 0.007261 

0 

0 


0.0566 


■ 

C 

6 . 44396 X 10" 32 

. .007601 

- 2 . 79674 X 10* 31 

-.032989 

. 165 X 10' 20 

.0103 


.0176 


Wm 

E 

0 

0 

- 4.75036 

-.056033 

0 

0 

-.279 

-.0174 

2 . 11667 X 10' 10 

4.0 

A 

0 

0 

- 2 . 06773 X 10' 32 

- 0.002439 

0 

0 

0.594 

0.0371 



C 

3 . 05200 X 10" 32 

-.003600 

- 1 . 41019 X 10' 31 

-.016634 

. 069 X 10"? 0 

.0043 

.179 

.0112 



E 

0 

0 

- 2.43143 

-.028680 

0 

0 

-.130 

-.0081 


a See fig. 1. 
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Component of electric dipole moment, or atomic units 








1 2 3 4 

I ntermolecular distance, R r , bohr 


Figure 5. - Derivatives of electric dipole moment with respect to molecular inter- 
nuclear distance for three configurations of a hydrogen molecule and a hydro- 
gen atom. Coordinates are shown in figure 1(a). Derivatives not shown are 
zero. 


Comparison with Other Investigators 


The accuracy of the molecular integrals was checked against Magnasco and Musso 
(refs. 41 and 42) among others. Comparison of the values of energy, average interac- 
tion energy, dipole moment, and dipole moment derivative with values of other investiga- 
tors is made in the following subsections. 

Energy. - The energy of Hg - H or has been calculated by a great many investi- 
gators. However, many of them used inaccurate molecular integrals. Others used 
over-simplified models. The most convincing check of the present calculation was made 
with a special linear symmetric configuration of Shavitt, Stevens, Minn, and Karplus 
(ref. 27, table XI) . Their model was the same as in this report except that they may not 

have used or thogonalized orbitals. They used a of 2. 27149X10 1 ® reciprocal meters 

in D 

(1.202 1/bohr) and £ =( of 1. 99937x10 xu reciprocal meters (1.058 1/bohr). They 

3 C IQ 

got an energy of -7.0213x10” joule (-1.6106 hartrees), which is exactly what the method 
of this report gave. 

A comparison of interaction energies for the linear configurations of this report 
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(r » 7.41x10“** m (1.40 bohrs)) is made with other investigators in figure 2(a). The 
energies of Hirschfelder, Diamond, and Eyring (ref. 4) are higher than this report be- 
cause they used the same value of orbital exponent for all three Slater -type Is orbitals. 
The energies of Conroy and Bruner (ref. 26) are lower, presumably because they in- 
cluded electron correlation in their wavefunction in a more realistic manner. The ener- 
gies of Shavitt, Stevens, Minn, and Karplus (ref. 27) are lower than this report because 
they used 15 Slater -type orbitals rather than three. The energies of Edmiston and 
Krauss (ref. 29) are lower than this report presumably because their molecular -orbital 
wave function, constructed from 36 s- and p-type Gaussian orbitals, was closer to the 
true wave function. In addition, two semiempirical calculations have been included, 

\ 

although a comparison with ab initio calculations is not entirely fair. The Porter and 
Karplus potential (ref. 2) correlates exchange reaction rates and trajectories fairly well 
(refs. 36 and 38) and is presumably realistic at inter molecular distances less than 
1.75x10“*° meter (3.3 bohrs). The Mason and Hirschfelder potential (ref. 33) is the sum 
of a first-order perturbation energy and two dispersion terms and is presumably real- 
istic for intermolecular distances greater than 2.65x10“*® meter (5 bohrs). 

A comparison of interaction energies for the scalene configurations is made in fig- 
ure 2(b). The Porter and Karplus potential (ref. 2) is again lower than the interaction 
energy calculated in this report. 

A comparison of interaction energies for the isosceles configurations of this report 
(r « 7.41x10“** m (1.40 bohrs)) is made with other investigators in figure 2(c) 
Hirschfelder 's energy (ref. 6) is higher than this report because he assumed all Slater- 
type Is orbitals to have orbital exponents of 1.88976x10*° reciprocal meters 
(1.0 1/bohr). Conroy and Bruner's energies (ref. 26) are lower, presumably because 
they included electron correlation in their wavefunction in a more realistic manner. 
Aroeste and Jameson's energies (ref. 15) are lower than this report, presumably be- 
cause they approximated some molecular integrals. Trivedi's results (ref. 31) are not 
shown because of only one-place accuracy for some of his three -center molecular inte- 
grals and other errors. The semiempirical potentials of Porter and Karplus (ref. 2) 
and Mason and Hirschfelder (ref. 33) are also shown. 

Average interaction energy from this report is compared with other investigators in 
figure 3. The semiempirical potential of Porter and Karplus (ref. 2) was averaged over 
H 2 orientations by use of equation (35) and is lower than this report because their poten- 
tial was lower for each configuration. The semiempirical potential of Mason and 
Hirschfelder (ref. 33) was averaged by means of the approximation 



int,A 


+ 



int,E 


(44) 
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which is fully analogous to equation (35). The interaction energy obtained by Amdur 
(ref. 34) from scattering of an H atom beam in Hg gas is the lowest of all for small in- 
termolecular distances. 

In summary, the interaction energy in this report is within the scatter of other ab 
initio calculations and has the proper dependence on inter molecular distance. At small 
intermolecular distances the dependence on H 2 orientation agrees with the Porter and 
Karplus potential (ref. 2). (It does not contradict Conroy and Bruner's finding (ref. 26) 
that the path of minimum energy involves a linear symmetric configuration because we 
fixed the H 2 internuclear distance at 7. 41599xl0”^ 1 m (1.401446 bohrs), but their saddle- 
point in their minimum energy path occurs for H 9 internuclear distance of 9.26047 
xlO m (1.75 bohrs).) At large intermolecular distances the dependence of interaction 
energy on H 2 orientation agrees with the Mason and Hirschfelder potential (ref. 33). In 
any case, the interaction energy is a sensitive function of the model used because it is a 
small difference between two large numbers (the first number being energy and the 
second being the sum of isolated H 2 and isolated H energies) . 

Dipole moment and its derivative . - There are no previous calculations of either of 
these quantities themselves, but for large enough intermolecular distance accurate re- 
sults may be obtained by classical methods: the H 2 quadrupole moment induces a dipole 
in the H atom. To draw meaningful conclusions from such a comparison, the compo- 
nents of the dipole moment should be expanded in normalized Legendre polynomials or 
normalized associated Legendre functions as follows: 

M z = Yj C z (R r )0^(cos 9) (45) 

l 

2- I/2 (M x * t4 y ) -££ Gj m <R r > ( W' :os e > eim4, < 46 > 

l m 

where 6 and 4> are shown in figure 6 . 


C 2 " 3 ' 


G 21 = -(6/5) 1/2 a Q/rJ 


In atomic units the classical method gives 


C 0 = ° 




(2/5) 1/2 a Q/ft* 


( 47 ) 


c 4 = o 
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Expansion coefficient, Cj or G21, 





where a is the polarizability of the H atom and Q is the scalar quadrupole moment of 
H 2 The value of a has been given by Pauling and Wilson (ref. 52), and the value of Q 
by Kolos and Wolniewicz (ref. 53), who use a different definition so that their values 
must be divided by two. The results are plotted for R r of 2. 38X10" 1 ® to 2..65X10" 1 ® 
meter (4. 5 to 5 bohrs) in figure 7. For comparison, the valence-bond results in figure 4 
were expanded by equations (45) and (46) and plotted in figure 7 for R of 0. 53X10" 10 
to 2. 12xi0" A meter (1 to 4 bohrs) . The valance -bond Cq and correctly approach 
0 as R r becomes large. The valence -bond C 2 is in good agreement with the classi- 
cal C 2 - However, there is a minor difficulty in the case of G 2 j: at large R r the 
valence-bond and classical values differ in sign, although both are small. This is be- 
cause the Is orbital centered on H in the valance-bond model cannot be polarized in the 
x direction. At smaller R r x-polarization of H will be approximated by changes in 
orbital populations. 

The derivatives of the components of dipole moment with respect to H 2 internuclear 
distance can be examined the same way. First they are expanded 


— = Y D,(R )© (cos 0) 
3r i 


(48) 


_3_ 

3r 


2' 1/2 <4 x + i4 y ) 


l m 


im<3> 


(49) 


In atomic units the classical method gives 


A 


D 0 = ° 


D 2 = 3 • (2/5) 1/2 aQ’/R^ 


d 4 = o 


J 21 = -(6/5) 1/2 aQ'/R* 


(50) 


where Q' is the derivative of Q with respect to H 2 internuclear distance. The results 
are plotted for R r of 2.38x10“^® to 2.65xl0 - ^ meter (4. 5 to 5 bohrs) in figure 8. For 
comparison, the valence -bond results in figure 5 were expanded and plotted in figure 8 
for R f of 0. 53X10" 1 ® to 2. 12X10" 1 ® meter (1 to 4 bohrs). The valence-bond Dq 
and correctly approach 0 for large. R r> The valence-bond D 2 is in excellent 
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1 2 3 4 5 


I ntermolecular distance, R r , bohr 

Figure 8. - Expansion coefficients of derivatives of components of dipole moment for H 2 -H. Valence 
bond values are plotted from 1 to 4 bohr. Quadruple induced values are plotted from 4. 5 to 
5 bohr. Dashed lines connecting corresponding coefficients were added from 4to 45 bohr for 
clarity. 

agreement with the classical D 2< However, for the reasons mentioned in the previous 
paragraph, the valence -bond J 21 and classical J 21 have different signs at large R r . 
Fortunately they are both small. 

CONCLUDING REMARKS 

Full valence-bond ab initio calculations of molecular hydrogen-atomic hydrogen 
(H 2 ~H) interaction energy and electric dipole moment were carried out with accurate 
values for all molecular integrals and optimization of each orbital exponent. Linear, 
scalene, and isosceles configurations were included with intermolecular distances from 
0.53x10"^ to 2.12xl0 - ^ meter (1 to 4 bohrs). A weighted average of the interaction 
energies of the three equilibrium configurations was then taken. 

The interaction energies obtained were within the scatter of other ab initio calcula- 
tions and had the proper qualitative behavior. However, they tended to be higher than 
other published values. 

Neither the dipole moment nor its derivative with respect to H 2 internuclear distance 
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had been previously calculated. The values obtained by the valence-bond model appeared 
to have the correct trends at large inter molecular distances and to be consistent with 
values obtained classically with one minor exception: the sign of the small transverse 
component was different for the valence -bond and classical methods. This was attribu- 
ted to a lack of polarizability of the H atom in the valence -bond model for large inter - 
molecular distances. 

The electric dipole moment and its derivative can be used to calculate pressure- 
induced vibrational absorption coefficients. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 8, 1973, 

503-04. 
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APPENDIX - SYMBOLS 


(abc) 

C 

Cl 

C, 


c pj 


ppj 


pq] 


D, 


,- 1/2 


< /2 


at 


E. 


E. 


int,A 


E 


int 


E 

m 


G ijki 


J l m 


mnop 


shorthand notation for antisymmetrized product of orthogonalized spin- 
orbitals 

column vector of coefficients of i b. 

] 

column vector of coefficients of cp. 
expansion coefficient for g z 

j th Cartesian coordinate of proton p 

1L 

j Cartesian coordinate of center of atomic charge p 
Cartesian coordinate of center of overlap charge pq 
element of C' 

expansion coefficient for 3u /Sr 
matrix of overlap eigenvalues 
inverse square root of d 
diagonal element of d - *^ 

electronic and nuclear repulsion energy of complex 
energy of atom 

electronic energy 

interaction energy for configuration A 
average interaction energy 

electronic and nuclear repulsion energy of diatomic molecule 
one-electron integral of orthogonalized orbitals 
expansion coefficient for E^ 
two-electron integral of orthogonalized orbitals 
expansion coefficient for 2 _1 ^(p x + i 
two-electron integral of Is orbitals 



H 

H 


H 


nk 



i 


J 


Zm 


K, 


kZ 


L. 


M 


jkp 

- 1/2 



m 


mane, 

rs 

m/3n/3 

Vs 


P’ 

mn 




Q 

Q' 



r 


Hamiltonian matrix with basis functions if/. 
Hamiltonian operator for electrons 
element of H 

Hamiltonian matrix with basis functions <p . 


element of H 


<P 


Planck constant divided by 2v 


(- 1 ) 


1/2 


expansion coefficient for d/dr |2“^^(|u i + in ) 

i x y 


kinetic energy integral 
nuclear attraction integral 
Lowdin transformation matrix 


element of M 


- 1/2 

X 


integer in expansions in normalized associated Legendre functions 
coefficients of spin-orbital products with two a spins 
coefficients of spin-orbital products with two (3 spins 
electron population in the orbital density X' m X' n irrespective of spin state 
electron population in the orbital density X D X a irrespective of spin state 

r , *x 

matrix of 

pq 

normalized one -electron population 

scalar quadrupole moment 
3Q/3r 

(Sp - ^q) R pq/ 2 (subscripts vary) 
distance between points p and q 

intermolecular distance (see fig. 1) 

internuclear distance of diatomic molecule 
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r ij 

distance between electron or proton i and electron j 

s 

matrix of 


pq 

V 

overlap integral for x n and x n 

u 

unitary matrix 

rs# 

U 

transpose of U 

V 

matrix relating and cp 

V 

transpose of V 

V i 

volume in configuration space for electron i 

w.,w.. 

expansion coefficient for 

W pq 

<?p + VV 2 

x,y,z 

Cartesian coordinates 

a 

polarizability 

a( i) 

spin eigenfunction of electron i with component of spin angular momentum 
along the axis of quantization equal to fi/2 

/3(i) 

spin eigenfunction of electron i with component of spin angular momentum 
along the axis of quantization equal to -K/2 

A pq 

displacement of overlap charge center from p-q geometric center, mea- 
sured in direction from p towards q and expressed in units of R 

r4 

? 

orbital exponent 

? P 

orbital exponent of x D 

r 

^po 

approximate optimum orbital exponent of x n 

0. 

l 

normalized Legendre polynomial 

0.. 

i] 

normalized associated Legendre function 

e 

angle between internuclear axis of molecule and line through atom bisecting 
the molecule 

4 i 

1L 

j n component of electric dipole moment 

r 

volume in configuration and spin space of three electrons 

$ 

azimuthal angle about z axis 
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<p 

Vi 

X 

xl 

x p a) 

X’ p (i) 

* 

± 

S2 


V? 


row vector with elements < cp^ 

antisymmetrized product of orthogonalized spin orbitals 

row vector with elements x p 
row vector with elements x' n 
12 orbital of electron i centered on proton p 

orthogonalized orbital of electron i 

wave function of system 

row vector with elements i/a 

formal orthogonalized valence -bond structure 

solid angle for molecule orientation 

Laplacian operator for the coordinates of electron i 
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